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The possibility of formation of a droplet phase (DP) inside a star and its consequences on the 
structural properties of the star are investigated. For nuclear matter (NM), an equation of state 
. (EOS) based on finite range, momentum and density dependent (FRMDD) interaction, and which 

' predicts that neutron matter undergoes ferromagnetic transition at densities realisable inside the 

, neutron star is employed. An EOS for quark matter (QM) with density dependent quark masses, the 

so-called effective mass model, is constructed by correctly treating the quark chemical potentials, 
i— H ' A comparative study of hybrid star properties as obtained within the usual bag model and the 

effective mass model shows that both these models yield similar results. Then, the effect of spin 
polarisation on the formation of DP is investigated. Using the EOS based on FRMDD interaction 
alongwith the usual bag model, it is also found that a droplet phase consisting of strange quark 
matter and unpolarised nuclear matter sandwiched between a core of polarised nuclear matter and 
a crust containing unpolarised nuclear matter exists. Moreover, one could, in principle, explain the 
mass and surface magnetic field satisfactorily, and as well allow, due to the presence of a droplet 
■ phase, the direct URCA process to happen. 
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!— h ' I. INTRODUCTION 

Since the discovery of pulsars (lj and their identification (5) as rotating neutron stars during the late 1960's, much 
theoretical work have been done to study the structure of neutron stars. At the core of these objects, the density 
could be as high as five to ten times the nuclear matter density (p ~ 0.16 fm~ 3 ) which falls down by several orders 
of magnitude in the crust region. The nature of matter at such high densities in the core is still an unresolved issue 
Over the two decades, many of the studies were addressed towards a resolution of this problem (see for e.g. 
It has been conjectured in particular that there exists a phase transition between nuclear matter (NM) and quark 
matter (QM) at high densities [5-11]. Wittcn [Q has further speculated that the strange quark matter (SQM) 
consisting of u, d and s quarks may be the absolute ground state of hadronic matter. If this is true, then the 
possibility that pulsars are objects made up of purely SQM cannot be ruled out. Recently, from a semiempirical 
analysis of the mass-radius relationship Li et al [ p"3] have suggested that Her X-l pulsar maybe a strange quark 
star (see also On the basis of these hypotheses, several authors [3,5-11] have studied the properties of stars 

consisting of only SQM, i.e. strange quark stars (SQS) and neutron stars having quark cores with nucleon envelopes, 
viz. hybrid stars (HS). In general, the transition between the nuclear and quark matter is considered to occur at 
a unique pressure determined by Gibbs criteria, ( we consider the transition to be first order). Inside the star, the 
pressure decreases smoothly outwards. The mixed phase has constant pressure independent of the volume fraction of 
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joint on the density profile of the star. Consequently, 
lq] has recently proposed that for systems with more 



the phases and therefore, the mixed phase collapses to a single 
the mixed phase cannot exist in the hybrid star. Glendcnning 
than one conserved charge, the QM and NM could coexist for a finite range of pressures. A well-defined mixed phase 
may therefore be present 15-17]] over a finite region within the star (hereafter referred to as the droplet phase ) 



Occurrence of a droplet phase (DP) have quite important observational consequences. In particular, the values of 
(3— equilibriated proton fraction x v in DP can be large enough to allow direct URCA process to happen, and hence, is 
of significance to the cooling process of neutron stars. However, whether or not such a DP is energetically favourable 
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depends largely upon the nature of dense matter. As yet, the exact nature of dense nuclear matter is not known, 
and several possibilities are equally probable. Very recently, we had shown jl8), using a well-defined, finite range, 
momentum and density dependent interaction, that pure neutron matter undergoes a ferromagnetic phase transition 
at density p ~ 4po- Using such an equation of state, we then made a detailed study of structural properties of neutron 
and hybrid stars, and showed that one could, in principle, explain both the surface magnetic field and the standard 
structural properties. An interesting result is that (3— cquilibriated proton fraction x p decreases sharply once the 
star matter becomes spin polarised (see Fig. 5 of Ref. One then naturally wonders, what is the effect of such 

phase transitions on the occurrence and properties of DP ? In particular, whether or not ferromagnetic domains, i.e. 
droplets of spin polarised matter, exist inside a star. And, what is the role of spin polarisation in the determination 
of x p , and thereby, on the possibility of direct URCA process. In view of these, we mainly investigate the following. 
Firstly, whether the proposed equation of state, based on an interaction with firm basis in the well known properties 
of nuclear matter and finite nuclei, which predicts a ferromagnetic phase transition permits a droplet phase. Secondly, 
if such a phase is allowed, what are the effects of spin polarisation on the star structure ? Finally, can one consistently 
describe the standard structural properties and surface magnetic field, as well as find the finite probability for the 
direct URCA process to happen ? 

Now, in addition to the properties of dense nuclear matter, the nature of the quark matter equation of state can 
also affect the possibility of formation of a droplet phase. The theoretical framework that has been used in general to 
study HS and SQS is the MIT bag model JlS[] . In this model, the non-perturbative aspects of quark-quark interaction, 
i.e. the confinement property is invoked through a bag parameter. An alternate way of implementing confinement 
has also been explored p0|-p2]|, which is hereafter referred to as effective mass model of confinement. In this approach, 
the quark masses are assumed to be density dependent. In the limit of high densities, their values tend towards the 



current quark masses; in the limit of zero density, the quark masses become infinite. Attempts have been made [£3 24 
to study the properties of SQS in this approach. In Ref. p3[ , it was noted that the properties of SQM obtained using 
the density dependent quark masses are quite different from those obtained using the MIT bag model. This was 



subsequently contested in Ref. |24|, where it was found that both the approaches yield similar results. An intriguing 
aspect in the latter study is that the pressure does not vanish at the density where the energy per baryon has a 
minimum. This contradicts the very definition of the ground state of a system. As shown in the following section, 
both these studies have incorrectly determined the quark chemical potentials. Due to the explicit density dependence 
of quark masses, there is an extra term, the so-called rearrangement energy, in the definition of chemical potential. 
Taking this term into account, it is shown that the energy minimum and zero pressure occur at the same density 
point. Further, we make a comparative study of properties of hybrid stars as obtained using the effective mass model 
and bag model of quark confinement. 

It may also be said here that in addition to the bulk properties of the equations of state, the finite size effects 
|jT6| , [25f such as Coulomb, surface and curvature contributions can affect the possible formation of the droplet phase 
inside a star. However, it is known that surface tension of QM is poorly determined Jig ], and is also found to be 
dependent upon the shape of the confining potential p6| . Moreover, the internal structure of the droplet phase is 
generally determined by the competition between the Coulomb and surface energies, and can adopt shapes like rods, 
plates and spheres p6| . Here, since we are mainly concerned with the effect of spin polarisation in the formation of 
the droplet phase, the finite size effects are presently not considered. 

II. THEORETICAL FRAMEWORK FOR QUARK MATTER EOS 

In the following, we briefly outline the procedure to obtain the quark equation of state in the effective mass 
approximation. 



A. The effective mass approximation 



The two important properties of QCD, namely the asymptotic freedom and the confinement, are adequately rep- 
resented by considering that quarks are bound by a local Dirac-scalar potential, S(r) ~ c/| r | J^ , ||[ . Presence of a 
strong Dirac-scalar component is found to be consistent with the concept of chiral symmetry p9[ . In terms of S(r), 
one can then define the quark effective mass as m q {r) = m q + c{\ r | p8pc| ]. Its average value in the bulk matter 
limit is taken to be of the form, 

m q =<m q (f)>=m° q + Cp^ n ; n = 1/3 > 0, (I) 
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where mj is the current quark mass, pb is the baryon number density and C and n are parameters to be determined 
by requiring that strange quark matter is stable or unstable against the normal nuclear matter. In this effective mass 
model, quark-quark interactions are assumed to be included through the density dependence of the quark masses. 
Further, Eq. ([lj) implies that in the asymptotically free regime, i.e. at high densities, m q — ► m q) the current quark 
mass. In the limit p& — ► 0, ra q — > 00; this gives rise to absolute confinement. 

B. Rearrangement energy and Hugenholtz-Van Hove theorem 

Consider a system of one kind of quarks with an effective mass m q defined in Eq. (Q). Then, the total quark number 
N q and the energy E q of the given system are determined using the expressions: 

JV,= J>(k), 

k 

E q = ^n(k)y/k? + m*, (2) 

k 

where k is the momentum and m q is in units of he, with K = c = 1. The occupation probability function n(k), at a 
given temperature T is 



n{k) = 



1 4- exp 



£fc - Pq 
T 



(3) 



where p q is the quark chemical potential. The single particle energies are determined using the standard definition, 
eft = dE q / dn{k). In the limit of zero temperature, one gets, 

e k = ^k 2 + m 2 + JL ^l[ mq k F ^k 2 F +m q - 
m„ In 



kF + \/kl + m 2 q 



= ^k 2 + m 2 q + U(p q ), (4) 

where p q (= 3pb) is the quark number density and the spin-colour degeneracy factor g = 6. The term U(p q ) in Eq. (^) 
arises solely due to the density dependence of quark mass. In the context of nuclear physics, such a term is normally 
referred to as the rearrangement energy |H],[32| ■ Because of this term, the single particle levels depend explicitly upon 
the value of the Fermi momentum kp. In the limit T — > 0, one can obtain from Eq. (Q) the chemical potential p q , 
equivalently the Fermi energy ep to be 



£f = ^ q |t=o= y kp + m q + U(p q ). (5) 
Similarly, the quark number density p q — N q /V and the energy density e q = E q /V are determined using Eq. (S) as 



Pq 



gkp 

6tt 2 ' _ 
, - [(2fc| + m 2 q k F )^k 2 F + m 2 




(6) 

where V is the total volume of the system. In the limit C — > 0, U(p q ) — ► 0, and one arrives at the usual Fermi 
gas model expressions for p q and e q . Starting with the standard definition (g) for energy, one thus arrives at a new 
feature in the expression for p q in the effective mass approximation as compared to the bag model. 

We now demonstrate that Eqs. (||) and (Jfy are consistent with the well-known Hugenholtz-Van Hove (HVH) 
theorem pS] arrived at in the context of interacting Fermi systems. The HVH theorem in general deals with the 
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single particle properties of an interacting Fermi gas at T = 0. For a system of one kind of particles, the theorem 
states that 



(ml 

It is straightforward to see that {dE q /dN q ) v — E q /N q + p q [d(E q /N q )/dp q } v . Then at equilibrium, i.e. at the density 
corresponding to zero pressure one has as a special case, p q — E q /N q . Thus, it relates the quark chemical potential 
at equilibrium to the average energy per particle. 



Using Eqs. (|) and (fjj), we then obtain, 

Pq 



k 2 F + m q + U(p q ), (8) 

which is same as in Eq. (|^) demonstrating the consistency with the HVH theorem. The results are e qually true for 
multicomponent quark systems also. We do this exercise in order to emphasize that earlier studies [ p3p4[| do not 
include such a rearrangement term in their definitions of p q . 



(dE, 
\dN q 



C. The equation of state 



Here, we construct, within the effective mass approximation, the equation of state of /3— equilibrated, electrically 
neutral quark matter. 

The total kinetic energy density of a system of non-interacting, relativistic quarks of flavour r and effective mass 
m T is given as, 



3 (m T c 2 ) 4 
8tt 2 (Tic) 3 



c rV /l + x 2 T (1 + 2x 2 T ) - ln{x T + y/l + xl) 



(9) 



where x T = p r F j (m T c), p T F is the Fermi momentum and is related to the quark number density p T of a given flavour as 

1/3 

p F = h(ir 2 p T ) . The densities pertaining to the three flavours can be expressed in terms of the total quark number 
density p q and the asymmetry parameters 5 ud and 8 US as: 



Pu = {Pq/3) [1 - <>ud - Sus] = (Pq/3)fu, 

pd = (p g /3) [1 + 2S ud - 5 US ] = (p q /3)fd, 

Ps = (Pg/3) [1 - $ud + 28 us ] = (p q /3)f s , 



(10) 



where 8 ud = (p d - p u )/p q , 5 US = (p s - p u )/p q and p q — p u + p d + p s = 3p b . Similarly, the energy density e L pertaining 
to a system of relativistic non-interacting electron gas is obtained as 



1 (m e c 2 ) 



8tt 2 (he) 



,^l + x 2 e (1 + 2a; 2 ) - ln(x e + ^l + x 2 e ) 



(11) 



where x e = p^ ie /(m e c), m e = 0.511 MeV, and pp. e is the electron Fermi momentum related to the density p e as 
p F , e = h(3iT 2 p e ) 1/3 . 

The equilibrium composition of the quark matter is then determined subject to the (3— equilibrium conditions, 



Mrf - Pu = Pe and p, d = Ps 



(12) 



and the charge neutrality condition, 



Pe = g (2p u - Pd - Ps) 



(13) 
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Using Eq. (|1C|) one obtains, p e — —(p q /3) (S u d + S us ). Similarly, one can express the chemical potentials p u , pd and 
p s in terms of the three quantities p q , S u d and S us as follows: 



Pu 



Pd 



Pie s 


\ Pic 


1 4- A , Pic 
L T u d OE q 


1 4- A Pic 
1 ~r Ous Ot q 






Pq dS ud 


Pq dS us 


<v 


) - 


1 - 5 ud de q 


S us de q 


dpdy 




Pq d5 ud 


Pq dSus ' 




1 - Oh. 


8 u d de q | 1 


- S us de q 


dp s/ 


j p„p d d p« 


Pq dS ud 


Pq d5 us ' 



\OPsJ PutPd OPq Pq OOud Pq Od U s 

where e g = e r is the total quark energy density ; r = u,d, s. The various energy derivatives occuring in Eq. ( |l4| ) 
are given by 
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d5 vd) Pq ,8 us 3 r T \ d5 ud 



where / T = p T /p b 



d5 ™)p q ,S ud 3 \ 5< W« 



t = 2^(^) 3 [avVT+^-M^r + + , (16) 



and X T = m r c 2 -\/l + a: 2 - 1S the chemical potential of the quarks with density independent masses which is also the 
usual bag model result. The chemical potentials corresponding to each flavour ( |l4| ) can also be expressed as follows: 

p T = \ T +Y / U T >. (17) 
t' 

Then using Eq. (|l4|), the (5— equilibrium conditions ( p"2| ) can be expressed in terms of the asymmetry parameters as 

9e " de J-=Q 



dSud d5 us 

2 de g 1 de q 

~ 'ax 1 ox - ~~ ( L °) 

Pq OOud Pq OOus 



where p e = <Jpp e (? + m^c 4: . ^From Eq. , the same can be reexpressed as 

Ad - X s — and X d - X u = Pe- (19) 

Thus, for a given baryon density pb = p q /3, the three quantities p e , 5 u d and S us are fixed by the Eqs. @ and @. 
The total energy density of quark matter cqm and the pressure Pqm for a given pb are then given as 

fQM = ^ e r(Pq, Sud, Sus) + e L(Pq, S u d, S U s), 

T 

-Pqm = P 9 CQM- (20) 

Opq 

Eq. ( pp| ) defines the equation of state of charge neutral, /3— equilibrated quark matter. 

In order to make a comparative study, we also consider the bag model picture. The quark matter equation of state 
within the bag model can be determined from Eq. (^TJ) itself by setting C = 0, and adding a bag energy density B to 
the total energy density cqm- One then gets, 

£ ba g = £qm[C = 0] + B, 

D _ £ bag /oi \ 

-Tbag — Pq -q Cbag- 

OPq 

Before, we discuss the results obtained within these two models, we give below the equations of state considered for 
nuclear matter. 
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III. NUCLEAR MATTER EQUATION OF STATE 



A. A simple parametrisation 

For the comparative study of properties of hybrid stars within the bag and effective mass models, we choose a simple 
parametrisation of nuclear matter EOS as given in Ref . [fL6| . The total energy density e„ for a system of neutrons 
and protons is taken to be of the form, 

e n = Pb [M n + E comp + E sym ] , (22) 

where M n is the nucleon mass. The compressional energy is given by -Ecomp = (Ko/18){pb/ Po — I) 2 and the symmetry 
energy is E aym — S (pb/ Po) (1 — 2a; p ) 2 , where x p is the proton fraction. Values for the nuclear matter incompressibility 
K , symmetry energy coefficient S and the normal nuclear matter density p are taken to be 250 MeV, 30 MeV and 
0.16/m -3 respectively [ p^[ . 

The optimum value of x p for a given pb is determined by the (3— equilibrium condition and the charge neutrality 
condition: 

Hn = P P + Pe and p e = Xpp b , (23) 

where p n and p p are the neutron and proton chemical potentials respectively. The total energy density enm and 
pressure Pnm corresponding to nuclear matter for a given baryon density pb is then given as, 

£NM = e n (Pb,X p ) + €L(Pb,X p ), 

v <9e N M , OA s 

-TNM — Pb~ £NM- 1^4 ) 

Opb 

The above equation defines the equation of state of the (3— equilibrated, charge neutral nuclear matter. 



B. Based on a realistic interaction with explicit spin degrees of freedom 

The phenomenological momentum and density dependent finite r ang e interaction employed here to obtain the 
equation of state is a modified version of Seyler-Blanchard interaction |32j| . To treat spin-polarised isospin asymmetric 
nuclear matter, the interaction has been generalised to include explicitly the spin-isospin dependent channel. The 
interaction between two nucleons with separation r and relative momentum p is given by, 



v eff(r,P,p) 



b 2 



d 2 {pi(ri) +P2(r 2 ))" 



-r J a 



r /a 



(25) 



where a is the range and b defines the strength of the repulsion in the momentum dependence of the interaction. The 
parameters d and 7 are measures of the strength of the density dependence, and p\ and p 2 are the densities at the 
sites of the two interacting nucleons. The subscripts r and s in the strength parameter C TS refer to the likeness I and 
the unlikeness u in the isotopic spin and spin of the two nucleons respectively; for example, Cu refers to interactions 
between two neutrons or protons with parallel spins, Ci u refers to that between neutrons or protons with opposite 
spin etc. The energy per nucleon E/A and the pressure P in the mean- field approximation can then be worked out 
IB3 as 



E/A=-Y J Prs 



' T h/2{Vrs) ^ 



1, 



m * v 1 ) + ~v° 

V TS y TSJ 1 2 



(26) 



TS 

Here, Jk{v) are the Fermi integrals, V® s and V 2 S are the single-particle and the rearrangement potentials, V^ s is the 
coefficient of the quadratic momentum dependent term in the potential and defines the effective mass m* s , T the 
temperature, and 77 is the fugacity given by r/ TS = (p TS — V® s — V 2 S )/T. For the unpolarised nuclear matter (UNM), 



2 ^3/2(777-3) 

3 Jl/iiVrs) 



+ v T ° s + ^(i-d\2 P y)v T \ + v?. 



(27) 



G 



the expressions for V® s etc are given in ref. |32fl. It is straightforward to extend these to the case of polarised nuclear 
matter (PNM) and are given in detail in Ref. |18| . 

One usually defines the neutron and proton spin excess parameters ( spin asymmetry ) as 

On = (PnT - Pnl)/p , 

Ctp = (p P T - p P i)/p , (28) 

where 

P = Pn + Pp = (Pn? + Pnl) + (PpT + Ppi) > ( 29 ) 

is the number density. We then define the proton fraction as x — p p / p. It is related to the isospin asymmetry 
parameter X as, 

X = (l-2x) = ( Pn - pp)/p . (30) 

We also define the spin excess parameter as Y = a n + a p and the spin-isospin excess parameter as Z = a n — a p . One 
can then express the energy per nucleon E/Aoi the NM at zero temperature as 

E/A — Ey + E X X 2 + EyY 2 + E Z Z 2 , (31) 

where terms higher than those quadratic in X, Y and Z are neglected. Here Ey is the volume energy of the symmetric 
nuclear matter, taken as —16.1 MeV and Ex is the usual symmetry (isospin) energy, taken to be 33.4 MeV. The 
quantities Ey and Ez are the spin and the spin-isospin symmetry energies of the NM respectively. We take |34|,[55],[l8| 
Ex = 33.4 MeV, Ey = 15 MeV and Ez = 36.5 MeV. Values of these coefficients are uncertain to an extent. 

The generalised hydrodynamical model of Uberall (3(| , which gives (E z /E x ) 1/2 ^1.1 fixes E z for a given value 
of Ex- The value of Ey is relatively more uncertain, and it was found in our earlier study on neutron stars 
that standard neutron star properties are best explained with Ey = 15 MeV and hence we use here the same set of 
parameters. The values of the parameters C TS , a, b, d and 7 are then determined by reproducing Ey, Ex, Ey, Ez, 
the saturation density of normal nuclear matter (po = 0.1533 fm~ 3 ), the surface energy coefficient (as — 18.0 MeV 
), the energy dependence of the real part of the nucleon- nucleus optical potential and the breathing- mode energies 
|p7[ . The parameters of the interaction so determined are listed below: 

Cu = -305.2 MeV a = 0.625 fm 

Ciu = 902.2 MeV b = 927.5 McV/c 

C u i = 979.4 MeV d = 0.879 fm 3n/2 

C uu = 776.2 MeV 7 = 1/6. 

With the above value of the parameter 7, the incompressibility of symmetric nuclear matter is K = 240 MeV. 

It may be said here that the above interaction reproduces quite well the ground state binding energies, root mean 
square charge radii, charge distributions and giant monopole resonance energies for a host of even-even nuclei ranging 
from le O to very heavy systems. We have also seen that for symmetric nuclear matter, our results agree extremely well 
with those calculated in a variational approach by Friedmann and Pandharipande(FP) with Vu+TNI interaction 
in the density range \p$ < p < 2p . But, for unpolarised pure neutron matter, the energies calculated [p| with our 
interaction are somewhat higher compared to the FP energies, particularly at higher densities. However, our results 



L8[ are very similar to those obtained from the recent sophisticated calculation of Wiringa et al |39| with UV14+UVII 



interaction, and hence the present interaction can be extrapolated with some confidence to neutron matter at high 



densities. Then, using this interaction it was found |18| that neutron matter undergoes a ferromagnetic transition at 
a density p ~ 4po- Before, we investigate the effect of such a ferromagnetic transition on neutron star properties, 
certain general aspects of SQM as obtained within the effective mass model of confinement is discussed below. 



IV. FLAVOUR SYMMETRIC STRANGE QUARK MATTER 



Having clearly specified, in the above two sections, the various equations of state considered for quark matter and 
nuclear matter, we are now in a position to make a detailed study of hybrid stars. To do so, firstly, we would like to 
demonstrate the consistency of our present QM calculation done within the effective mass model, and as well make 
a preliminary study of the general aspects of the strange quark matter equation of state given by Eq. ( pp[ ) . For this 
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purpose, we take the current quark masses of u, d and s to be zero. As the parameters n and C (defined in Eq. ([[]) 
by the density dependence of the quark masses) are flavour independent, the quark masses are the same and hence, 
the matter is flavour symmetric, i.e. f u = f d = f s = I. 

The mass parameters n and C are chosen so that at the density where pressure P = 0, (M / A) p=0 < 930 MeV. Here 
M/A is the total mass per baryon of the flavour symmetric SQM. In Table 1, the values of (M/A) p _ Q so-obtained 
are shown for several values of n and C. The main result of this preliminary investigation is shown in Fig. 1, where 
we have plotted M/A and P as a function of the baryon density pb- Values of n and C are taken to be 2/3 and 
85 MeV jmT 2 . As expected, the minimum of M/A and P = occurs at the same density point for all values of n 
and C considered here. It may be noted that in the study in Ref. p4j ], the above mentioned density points did not 
coincide; this could be due to the noninclusion of the rearrangement term in the single particle energy. The role of 
this rearrangement energy in quark chemical potential is illustrated in Fig. 2, where the ratio pd/^d is displayed as 
a function of baryon density pb- The chemical potentials Xd and p,d are as defined in Eq. (fT^). The deviation from 
unity in the case of the effective mass model is due to the rearrangement term. The effect of this term is dominant 
in the low density regime. As the density pb increases, pd/Xd — ► 1, the bag model result; this is due to the fact that 
as pb — > oo, m T — ► m®, and hence the rearrangement term U T — ► as can be seen from Eq. (IHf). Similar curves 
for other flavours are not shown as p s /X s = pd/^d = Pu/X u . 

We have also studied the dependence of the properties of strange quark star on the mass parameters n and C. The 
total mass and the radius of the star are obtained by solving the general relativistic Tolman-Oppenheimer-Volkoff 
(TOV) equation @, 



dP{r) G [e(r) 


f P(r) 


[m{r)c 2 + 47rr 3 P(r)] 


dr c 4 




\ 2Gra(r)" 
-L — *> — 





(32) 



where, 

m{r)c 2 = f e{r')d 3 r'. (33) 
J o 

The quantities e(r) and P(r) are the energy density and pressure at a radial distance r from the centre of the star, 
and are given by the equation of state (pQ). The mass of the star contained within a distance r is given by m(r). 
The size of the star is determined by the boundary condition P(R) = and the total mass M of the star integrated 
upto the surface R is given by M — m(R). The single integration constant needed to solve the TOV equation is P c , 
the pressure at the center of the star calculated at a given central density p c . The values of the total mass M max 
and the radius R of the star corresponding to the maximum mass configuration are given in Table 1. The central 
density p c in units of p is also given. In order to understand the behaviour of the various structural properties of SQS 
with respect to the change in the mass parameters n and C, we have calculated the incompressibility K q — (dP/dp q ) 
of the star matter, which is displayed in Fig. 3 as a function of pb for three sets of (n, C). It is found that with 
increase in n (keeping C fixed), the EOS becomes stiffer. With n fixed, if C is increased, the EOS becomes effectively 
a little stiffer(averaged over the density of the star) though at low densities, it is somewhat softer. ^From Table 1, 
we note that as the EOS becomes stiffer, (M/A) p _ n and p c increases whereas R and M max decreases. These results 
corroborate the findings by Haensel et al [M in the bag model. We expect these features to be valid even when the 
/3— equilibrium conditions are considered, since it is in general found QlOfj that (3— equilibrated SQM is nearly flavour 
symmetric and the electron density is negligibly small. 

Finally, taking into account both 0— equilibrium and charge neutrality conditions given by equations (]lj) and ( |T^ ) 
respectively, we determined the mass-radius relationship of the strange quark star by solving the TOV equation in 
both the models described by Eqs. (|2(]) and ((H). Results so obtained are displayed in Fig. 4 and it can be seen that 
results obtained within the two models of quark confinement are quite similar. It may be stressed here that the values 
of -Mmax and R are dependent on the model parameters. It is found that with suitable choice of the parameters, 
it is possible to obtain the maximum mass and the radius of the star well within the acceptable limits. For the 
sake of comparison, we have shown the results obtained for a neutron star in the same figure. The structure of the 
neutron star is determined using a composite EOS, i.e. Feymann-Metropolis-Teller |ic[| , Baym-Pethick-Sutherland 
pl| , Baym-Bethe-Pethick [ fill a nd that given by Eqs. (|26[-p7|) with progressively increasing densities. It is already 
known from previous studies |lO| that curvature of the M — R curves in the low mass region as obtained from the bag 
model for SQS is opposite to that of the neutron star. This difference is also borne out in the present calculation using 
effective quark masses. Thus, our present calculation of quark matter EOS within effective mass model is consistent 
with the general findings in literature. 
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V. HYBRID STARS : A COMPARATIVE STUDY OF TWO MODELS OF QUARK CONFINEMENT 



Strange quark matter may not be the absolute ground state of hadronic matter. If this is so, then it is possible 
that neutron stars are not made up of only quark matter, but have a quark core enveloped by nuclear matter, usually 
referred to as hybrid stars. Since the quark-hadron transition is normally taken to be of first order, there exists a 
mixed phase consisting of droplets of quark matter and of nuclear matter in thermodynamic equilibrium. As this 
occurs at a unique pressure, such a phase cannot exist inside the hybrid star. Taking into account the recent viewpoint 
of Glendenning [lq| as already mentioned in the introduction, the existence of this droplet phase inside the hybrid 
star cannot however be ruled out. In this section, we investigate the structure of this droplet phase and its observable 
consequences on the structural properties of the neutron star using the two models of quark confinement, viz bag and 
effective mass models. We thereby perform calculations similar to the earlier studies J24|,|23|], but with the correct 
treatment of qua rk chemical potentials. Moreover, we make use of the simple parametrisation of nuclear matter EOS 
given by Eq. (|24| ) in order to be able to directly compare our results to those obtained in earlier studies [^6| , p3| . 

A. Equation of state for droplet phase 

We consider the droplet phase to be present in a uniform background |l(J of electron gas. The system as a whole 
is charge neutral, and each phase of matter is subjected to appropriate /3— equilibrium conditions. 

For the nuclear matter phase, we consider the EOS as given by Eq. (|24|). The appropriate 3— equilibrium condition 
is /j, n = Hp + /i e . The total charge density pf pertaining to the NM phase is given by p^ — p p — p e = x p p^ — p e , 
where p v is the proton density and p^ is the total baryon density in the NM phase. For given values of p^ and 
electron density p e , the proton fraction x p is determined from the 3— equilibrium condition. In the case of quark 
matter phase, the equations of state are given by Eqs. (^) and (|2l]). The total charge density p® in this phase is 
p® = (2p u — pd — Ps)/3 — p e , where the flavour densities are as defined in Eq. (|Io|). For given values of p?(— p q /3 = 
(Pu + Pd + Ps)/3) and p e , the asymmetry parameters 5 u d and 8 US are determined from the B— equilibrium conditions 
given by Eq. (|l8| ) . 

The two conserved charges of this uniform droplet phase are the total baryon number and the total charge. The 
corresponding baryon density pb and charge density p z are expressed in terms of the volume fraction \ as: 

p z = X P? + (1 - x)Pz, (34) 
Pb = XP? + (1 - x)Pb- (35) 

Global charge neutrality requires that p z = 0. Thus for given values of pb, x and Pb > the electron density p e is 
fixed by charge neutrality condition given by Eq. (H) and p? is determined from Eq. (|3|). To arrive at the EOS 
of the droplet phase at a given density pb, x an d p^are to be determined from mechanical and chemical equilibria. 
The chemical equilibrium between the NM and QM phases requires that p n = 2/id + Pu and p p = pd + 2p u . From 
mechanical equlibrium, one has Pqm = -Pnm, where the pressures in the two phases are as defined in Eqs. ( p4[ ) and 
@|!|). The total energy density £t of the droplet phase is given by 

(t = X e QM + (1 - x)cnm, (36) 

where 6qm and cnm are calculated with the constraints of mechanical, chemical and 3— equilibrium conditions as 
well as with the condition of global charge neutrality. 

For the study of the hybrid stars, the parameters in the effective mass model and the bag model are chosen so that 
(M/A) p=0 > 930 MeV. The parameters chosen in the effective mass model are n = 2/3, C — 125 MeV fm~ 2 and 
m° = 250 MeV. In the bag model, we have taken B 1 ^ = 175 MeV and m° = 250 MeV. 

B. Results and Discussions 

The droplet phase obtained with global charge neutrality condition (as detailed in Sect. V.A.) is found to be 
energetically more favourable than the mixed phase determined by the Gibbs criteria, where each of the two phases 
is separately charge neutral. This is illustrated, in the effective mass and bag models, in Figs. 5 and 6 respectively. 
The dashed line corresponds to the mixed phase as obtained by the common tangent method. It is seen that the 
droplet phase extends well beyond the mixed phase (dashed line). This was also noted in earlier studies Jl^,|l^]. In 







the limits of low density (\ — > 0) and high density (x — > 1), the droplet phase smoothly joins with the NM and 
QM equations of state respectively as it should. 

It may be parenthetically noted that the "down-turn" behaviour in the quark matter EOS in the low density region 
in Fig. 5 is due to the choice of the parameter n. In the limit pb — ► 0, the quark energy density e q [Eq. (P)] behaves 
as, 

5/3 

e 9 ~ h m q p q , 

m q 

~p n b +5/3 + C P l n . (37) 

If n — 1, in the limit p b — > 0, e q — > C in agreement with the usual bag model result. For n = 2/3, in the limit 
Pb — ► 0, e q — > 0, and this causes the "down-turn" in the quark matter EOS. We have chosen n = 2/3 as the 
Dirac-scalar potential [Eq. (0)] is generally taken to be of harmonic oscillator type, i.e. I = 2. 

We now explore the possible consequences of the presence of a droplet phase on the structure of the hybrid star. 
To do this, we solve the TOV equation using the EOS shown in Figs. 5 and 6. The baryon number density and 
mass distributions inside the star, corresponding to the maximum mass configuration, determined using the effective 
mass and bag models are shown in Figs. 7 and 8 respectively. The hybrid star structure in the effective mass is 
characterised by the central density p c = 0.83 /to -3 , the radius R = 11.8 km and maximum mass M max = 1.48Mq. 
In the bag model we have p c = 0.88 /m -3 , R = 11.6 km and M max = 1.52M Q . One sees from Figs. 7 and 8 that the 
droplet phase extends in both the models from ~ 1 km to ~ 8 km, and stars contain a core of about 1 km made of 
only QM and a crust region of about 4 km made up of onl y N M. We find that the droplet phase encompasses most of 
the volume of the star in agreement with an earlier study |ll| . In a recent study jl7| using the color-dielectric model, 
it was however noted that the substantial part of the core of the star contains only QM, and the droplet occupies 
roughly 40% of the total volume. It must be stressed that these finer details are in general model-dependent. Finally, 
the main conclusion of this section is that with the correct treatment of quark chemical potentials, we find that the 
hybrid star properties as obtained within the bag and effective mass models are quite similar. For this reason, for the 
study on the effect of spin polarisation on droplet phase, done in the following section, we use only the standard bag 
model. 

VI. SPIN POLARISED NUCLEAR MATTER AND HYBRID STARS 

In this section, we explore the properties of the hybrid stars using the proposed NM equation of state given by Eqs. 
( p6| ) and fl2^), which predicts a ferromagnetic phase transition at high densities. 

A. Spin polarisation and droplet phase 



Firstly, to know whether there exists a phase transition between the nuclear matter phase described by Eqs. 
and the QM phase characterised by the bag model EOS, Eq. (|2l]), we compare the total energies per baryon obtained 
in the two phases. The parameters chosen for this purpose are as given in section III.B for NM phase, and for QM, 
we have taken B 1 / 4 = 180 MeV and m s = 250 MeV. Further, in the case of NM, the total energy is minimised with 
respect to both the neutron and proton spin polarisation parameters, a„ and a p respectively. This energy minimised 
nuclear matter calculation is hereafter referred to as EMNM. In addition, we also calculated the total energies for 
unpolarised nuclear matter (UNM), where a n — and a p — at all densities. Results of these calculations are 
displayed in Fig. 9. 

One can observe a slight bend in the EMNM curve at p ~ 0.7 /to -3 . This point indicates the ferromagnetic phase 
transition density, pfm, where the optimum value of a n sharply increases from zero to unity. ( This feature is clearly 
illustrated in Fig. 3 of Ref. |18[ ). ) It is also seen that a polarised nuclear matter (PNM) phase is energetically 
favourable over the density region pfm < Pb < /°hq, where phq — lOpo is the NM-QM phase transition density. It 
would be interesting to know whether there exists a droplet phase consisting of droplets of SQM and PNM. For this 
purpose, we made an attempt to solve for the droplet phase following the same procedure detailed in section V. It was 
then found that there exists no solution for phase equilibrium between SQM and PNM phases. In order to understand 
this result, we do the following analysis. 

Firstly, to show the effect of increasing proton fraction on a„, we calculated, for given values of x p and pb, the 
optimum value of a n by energy minimisation. Values of a n so determined are plotted in Fig. 10 as a function of 
density for few values of x p . ( It may be said here that for x p < 0.5, values usually found in the droplet phase, the 
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maximum of the optimum value of a p , obtained with increase in density, is very nearly zero. ) As x p increases, the 
value of Pfm increases and the maximum value of a n , given by (1 — x p ), decreases. Therefore, in general, one might 
say that NM phase with large value of x p disfavours spin polarisation. On the otherhand, as it was noted |Q before, 
as NM gets spin polarised the 0— cquilibriated x p sharply falls to zero, which is also evident from Fig. 11, where we 
have shown mass per baryon versus density for certain cases. One can see that there is a slight deviation, for densities 
Pb < Pfm, between the EMNM (with x p — 0) curve and the EMNM curve obtained taking 0— equilibriation into 
account, which is essentially due to non-zero, though small, values of x p . However, as density is increased further, 
the two curves coincide exactly illustrating the fact that 0— equilibriated x p values become zero at p ~ pfm- Another 
interesting aspect is that beyond the ferromagnetic transition density, the PNM phase is energetically more favourable 
than the UNM with large values of x p . Hence, for densities pb > pfm, droplets of UNM with x p > would prefer to 
become PNM with x p = 0. That is, one expects a transition from a state (x p > 0, a n — 0) to a state {x p — 0, a n ~ 1). 
But, this change of state, in the context of formation of a DP, is disallowed for the following reasons. 

In contrast to the usual mixed phase defined by the common tangent method, in the case of a droplet phase, the 
electron density is same both in NM and QM phases. Restricting to densities pb > pfm, the following three cases are 
considered. 

CASE 1 : Firstly, let us consider the simplest situation with p e = and N e = 0, where we have neglected the small 
er mass. Then, the optimum configuartion for QM phase gives p u = Pd = Ps and N u = Nd = N s , where we have 
taken m u = rrid = m s . For the NM phase, at pb > pfm, the optimum configuartion gives x p = and ct n = 1. Thus, 
one finds that the two phases are separately charge neutral, and hence corresponds to the usual mixed phase with 

iV e = 0. 

CASE 2 : Secondly, let us consider as in case 1, p e = and N e — 0. But now, we insist that m s > — m u . Then 
the optimum configuration for the QM phase gives p u = pd = Ps with N s < Nd = N u , and thereby the QM droplets 
are now positively charged. For the global charge neutrality condition to be satisfied, the NM droplets should be 
negatively charged. However, according to the optimum configuration, the NM droplets consists of only spin polarised 
neutrons. Therefore, in this case, there exists no DP solution. 

CASE 3 : Finally, let us consider the general situation, where p e ^ ( and N e ^ ), and m s > nid — m u . Now 
with finite electron density, the 0— equilibrium condition pertaining to the NM droplets, at a density pb, is given by 
(1/ pb) dtn/dxp + p e = 0. Since /i e > 0, the slope factor de n /dx p must be negative. From Fig. 12, it can be seen 
that this is satisfied for the configuration UNM, i.e. x p ~ 0.5 and a n — 0. This then suggests that NM droplets 
with pb > Pfm i n DP may remain unpolarised with x p > 0, which however is not the optimum configuration for 
NM phase. Thus, due to the nature of droplet phase, that is, constancy of electron density in both phases, and the 
interplay of x p and a„, one finds that NM droplets with pb > pfm in DP may remain unpolarised with x p > 0. It 
may be remembered that the usual mixed phase between SQM and PNM however exists. In the following subsection, 
we explore the possibility of forming a DP with UNM phase and SQM phase. 

B. Properties of hybrid stars 

^From our earlier study on hybrid stars, it can be seen ( from Table 2 of Ref. |l^] ) that for the chosen set of 
parameters, the central density of the star is about 6po- Therefore, the UNM-SQM phase transition at p ~ 0.8 /m -3 
is of more relevance. Following the same procedure described in section V, we then solved the set of simultaneous 
equations corresponding to 0— equilibrium, charge neutrality, chemical and mechanical equilibria. Results obtained 
for the droplet phase along with the usual mixed phase are shown in Fig. 13. It can be seen that as in other cases, 
the droplet phase extends well beyond the mixed phase. Interestingly, the droplet phase is energetically favourable 
even than the EMNM case over the region 0.5 < p < 0.9. Further, in this particular case, only a part of DP is allowed 
as EMNM configuration takes over at high densities. Thus, one has according to the minimum energy criterion, a DP 
sandwiched between UNM and PNM phases. We now explore the possible consequences of the presence of a DP on 
the structure of hybrid stars. 

In Fig. 14, we show the proton fraction obtained in the DP u sing the realistic interaction [Eq. (|25|)], and compare 
with the one obtained within the simple parametrisation [Eq. It can be seen that the general trend is quite 

similar in both the cases. Values of x p as obtained with the realistic interaction is somewhat lower than the simple 
one. As the bag model parameters are more or less the same, this could be due to the fact that the EOS based on 
realistic interaction is relatively softer than the simple parametrised EOS. Nevertheless, values of x p is high enough 
to allow for direct URCA process to happen. We now solve the TOV equation (|32| ) using the EOS shown in Fig. 
13, and obtain the structural properties such as mass and size of the star. The baryon number density and mass 
distribution inside the star, corresponding to the maximum mass configuration, are plotted in Fig. 15. The maximum 
mass configuration is characterised by the central density p c = 1.02 /m -3 , the radius R = 11.4 km and maximum 
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mass M max = 1.62M©. Now, from Fig. 15, one observes that there exists a mixed phase consisting of droplets of SQM 
and UNM sandwiched between a PNM core and UNM crust. The droplet phase extends from ~ 2 km to ~ 6 km. 
More interestingly, there exists also a core of about 1 km consisting PNM. Thus, the salient features of this star 
structure are the following. Firstly, the mass and size of the star obtained in the present calculation is well within 
the acceptable limits. Secondly, due to the presence of a droplet phase, the value of x p found inside the star allows 
direct URCA process to happen, and thereby have important consequences on the cooling process of neutron stars. 
Finally, due to the presence of a spin polarised core, we could also, in principle, explain the surface magnetic field of 
neutron stars. Value of the magnetic field at the surface of the star due to the presence of a ferromagnetic core can be 
estimated using the relation H ~ a n p, n N/R 3 , where N ~ 1.4 x 10 54 is the number of neutrons in the ferromagnetic 
phase and R ~ 11.4 km is the star size. Then using | p n |= 1.91 nuclear magneton, we get H ~ 10 13 G. Thus, the 
surface magnetic field is somewhat overestimated compared to the observed value ~ 10 12 G. However, by fine-tuning 
the parameter Ey, as pfm is sensitive to this parameter, one could arrive at the right value. 

Before summarising our findings, we would like to stress that the star's internal structure very much depends upon 
the model parameters chosen. For example, keeping Ey fixed at 15 MeV, and decreasing the bag parameter from 180 
MeV to 170 MeV, one can see from Fig. 16 that PNM phase is no more allowed by energy minimum criterion. On 
the otherhand, with the bag parameter fixed at 180 MeV and increasing Ey from 15 MeV to 17 MeV, one finds that, 
though there exists a small density region where DP is favoured, PNM phase occupies larger phase space than the DP 
as illustrated in Fig. 17. Hence, in this case, the surface magnetic field will be largely overestimated. However, within 
an allowed set of parameters, we could, in principle, quite satisfactorily explain the structural properties, surface 
magnetic field as well as allow for the formation of a DP. 

VII. SUMMARY 

To summarize, firstly, we have constructed a quark matter equation of state in the effective mass approximation 
and then applied it to investigate some properties of strange quark matter with particular reference to hybrid and 
strange quark stars. In this effective mass model, the quark masses are chosen to be density dependent and the absolute 
confinement of quarks as well as the asymptotic freedom are adequately taken care of through this density dependence. 
With these effective masses, it is seen that the chemical potential acquires an extra term, the so-called rearrangement 
term. Presence of such a term in the single particle energy is consistent with the well-known Hugenholtz-Van Hove 
theorem. Thus, the present calculation, in contrast to earlier studies, correctly treats the quark chemical potentials 
in the effective mass model. 

Secondly, we make a comparative study of hybrid star properties as obtained within the effective mass and bag 
models. For nuclear matter equation of state, we first used a simple parametrisation. It is found that a well-defined 
mixed phase of droplets of quark and nuclear matter sandwiched between a pure quark matter core and a crust of 
only nuclear matter exists. This mixed phase is found to occupy most of the volume of the star. The main result of 
this comparative study is that both the bag and effective mass models yield similar results. 

Thirdly, using a realistic equation of state of nuclear matter based on a finite range, momentum and density depen- 
dent (FRMDD) interaction, which predicts a ferromagnetic phase transition at a density p\> ~ 4,00, we investigated the 
effect of spin polarisation on the formation of a droplet phase. It is found that as nuclear matter gets polarised, the 
value of proton fraction x p drops sharply to zero, or vice versa, as x p increases, the transition density p-pyv increases 
and hence matter with x p > prefers to remain unpolarised. Because of this interplay between a n and x p , there 
exists no mixed phase consisting of droplets of strange quark matter and polarised nuclear matter. 

Finally, using the same equation of state based on FRMDD interaction alongwith the usual bag model, we find that 
a droplet phase consisting of strange quark matter and unpolarised nuclear matter can however exist. This droplet 
phase is sandwiched between a core made up of polarised nuclear matter and a crust containing unpolarised nuclear 
matter. This, in principle, allows for a satisfactory explanation of the structural properties and surface magnetic field 
of neutron stars and keeps room for their rapid cooling due to the direct URCA process. 
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FIG. 1. Total mass per baryon M/A and total pressure P corresponding to flavour symmetric strange quark matter is 
shown as a function of baryon density pb in the effective mass model. M n is the nucleon mass and p is the normal nuclear 
matter density. 
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FIG. 2. Quark chemical potential [see Eq. (y_7|)] obtained in the effective mass and bag models are shown as a function of 
baryon density ph. 
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FIG. 3. Incompressibility K q of flavour symmetric strange quark matter is shown as a function of baryon density pb in the 
effective mass model. Values of the mass parameters n and C chosen are given. p is the normal nuclear matter density. 
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FIG. 4. Mass-radius relationships of strange quark stars (SQS) as obtained in effective mass and bag models are compared 
with a neutron star (NS) one. 
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FIG. 5. Total energy densities of nuclear matter (NM), quark matter (QM) and droplet phase are shown as a function 
of baryon density pt in the effective mass model. The dashed line represents the mixed phase obtained by common tangent 
method. M n is the nucleon mass. 
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FIG. 6. Same as in Fig. 5, but for the bag model. 
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FIG. 7. Baryon density [p(r)] and mass [m(r)] distributions corresponding to the maximum mass configuration are shown 
as a function of the radial distance r in the effective mass model. M sun (= Mq) is the solar mass and p c is the baryon density 
at the core (r = 0) of the star. 
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FIG. 8. Same as in Fig. 7, but for the bag model. 
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FIG. 9. Total mass energies per baryon of quark matter (QM), nuclear matter minimised with respect to spin polarisation 
parameters (EMNM), unpolarised nuclear matter (UNM) and mixed phase are shown as a function of baryon density pb- 
Parameters chosen are : B 1/4 = 180 MeV, m s = 250 MeV and E Y = 15 MeV. 
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FIG. 10. The value of spin polarisation parameter a„ obtained from minimisation of energy for particular values of x. 
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FIG. 11. Mass energies per baryon obtained for pure neutron matter (x p = 0) and /3— equilibriated nuclear matter (solid 
line) considering energy minimisation with respect to spin polarisation parameters are shown. Results obtained for unpolarised 
nuclear matter (UNM) with particular values of proton fraction x p are also plotted. 
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FIG. 12. Total energy density e n of nuclear matter obtained considering energy minimisation with respect to spin polarisation 
parameters are shown as a function of proton fraction x p for three particular values of baryon density pi,. 
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FIG. 13. Same as in Fig. 9, but the droplet phase consisting of unpolarised nuclear matter and strange quark matter 
droplets are also shown. 
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FIG. 14. Values of proton fraction a;„obtained for droplet phase using bag model equation of state (EOS) for quark matter, 
and with a simple parametrised EOS ( pi] ) and a realistic one ( p6|j2^ ) for nuclear matter are shown as a function of volume 
fraction \- 
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FIG. 15. Baryon density [p(r)] and mass [maldistributions corresponding to the maximum mass configuration obtained 
with bag model and realistic nuclear matter ( pilJi^ ) equations of state are shown as a function of the radial distance r. 
M sun (= M©) is the solar mass and p c is the baryon density at the core (r = 0) of the star. (For details see section VI.) 
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P b (fm 3 ) 

FIG. 16. Same as in Fig. 9, but with B 1/4 = 170 MeV. 
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TABLE I. Values of mass per baryon of flavour symmetric strange quark matter at zero pressure(M/y4) p=0 , and central 
density p c , radius R and mass M max corresponding to the maximum mass configuration of the strange quark star are given for 
few sets of the mass parameters (n, C). p is the normal nuclear matter density. 
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